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ARSTRACT 


Computations indicate that when utilized in an idealized 
ranjet at a flight Mach number of 2.5 and an altitude of 
37,500 feet, a constant area axial flow nuclear heat exchanger 
will be expected to produce 9,500 pounds of thrust with an air 
inlet Mach number of 0.2, a weighted radial temperature aver- 
age of 20009? and a frontal area of 7,7 square feet. Molyb- 
denum alloy with a disilicide coating was selected as the 
container material. Sandwiched between the molybdenum was a 
fissionable alloy of highly enriched uranium and zirconium. 
From a nuclear aspect the unreflected heat exchanger becomes 
critical at a diameter of 0 inches and a length of 60 inches 
with a totsl weight of 8,052 pounds which includes 3,622 
pounds of 90 per cent enriched uranium. 
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INTRODUCTION 


That a constant-area axial-flow compect nuclear heat 
exchanger as envisioned in Figure 1 when installed in an 
idealized ramjet would produce significant thrust was the 
basis of this thesis. Feasibility conclusions were limited 
to the information revealed by heat transfer and nuclear 
analysis under one assumed steady state condition. The 
problems of nuclear control, shielding, thermal stresses, 
and performance under different inlet Mach numbers were not 
considered. 

The nuclear heat exchanger as shown in Yiguro 2 incor- 
porates thin coneentrie cylinders which are the solid ris- 
81onable containers. The alloyed fizsionable phase is 
sandwiched between two layers of non«fissionable structural 
material so that the enriched uranium fuel cores are sur- 
rounded by a continuous network of structural metal for 
strength and containment at high temperatures. Air flows 
within the annuli and acts as the coolant and thrust arent. 
It is recommended that the nuclear flux be flattened by a 


eentral control core and a peripheral reflector; however, the 


beneficial effects of a reflector and central control core 
were not considered, 

From the assumed inlet condition, various maximum eyl- 
inder or wall temperatures were selected arbitrarily for 


thrust resulte. That the high temperature fuel elements and 


2 


structural materials were possible metallurgically were next 
investigated. Whether or not an unreflected homogeneous 
reactor of this configuration can be made eritical was deter- 


mined by using one group diffusion theory. 
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LIST OF SYMBOLS 


Cross section area normal to flow, (εὐ), or the mass 
number of an element (gm/rm atomic weight). 
Wetted area, (ft*). 
Buckling, (emf). 
Local velocity of sound, (ft/sec). 
Specific heat at constant pressure, (BTU/Lb Y), 
Diffusion coefficient fer flux, (em). 
Hydraulic diameter, D, - D,, (ft). 
Inside diameter, (ft). 
Outside diameter, (ft). 
Neutron leakage, (n/om? soc). 
Friction factor defined as 2g T. ov®, 
Acceleration of gravity (32.2 ft/sec^). 
Mass velocity (1b/see rt). 
Coefficient of heat transfer between fluid and surface, 
(BTU/hr ft? op), 
h, based on (Tow - Σο aw)? 
h, based on logarithm mean of LA - T). 
h based on T ow è T). 
h, based on (To ow * To)’ 
Extrapolated distance, (em). 
Neutron multiplication factor. 


Thermal conductivity, (BTU/hr rt? °F /ft). 
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L Length of heat exchanger heating element in the 
longitudinal direction, (ft). 

Mass rate of fÓ1ow, (1b/sec). 

Mach number or molecular weight (gms/gm molecular 


JR 


weight). 
Number of neutrons per em. 
Number of nuclei por om^. 
Avogadros number (6.02 x 102? nucleií/gm atomic weight). 
R Prandtl number, σομ/Κ. 


mm uw 5 


E Reynolds number, evD, /h. 
Stanton number, he / C90 


a το 


Absolute stream pressure (1b/1n*). 
3 Isentroplc stagnation pressure (1b/1n*). 
Heat-transfer rate (BTU/hr). 


6 ο vw 


Volumetric heat source which 1s equal to Eg,j, 
(BTU/sec em). 

Radius (ft). 

Absolute stream temperature (?R). 

Thrust (1b). 

Stagnation temperature (°R). 


Qu o" 


μα 


6, aw Adlabatic stagnation wall temperature (°R). 
Stagnation wall temperature (°R). 
Velocity (ft/sec). 

Volume (re? or en). 


Weight (1b). 
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Distance along heat exchanger heating elements (ft). 
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Subscripts 

b Bulk (T). 

f Film t, = (t + 4 )/2. 
W wall. 

0,1,2 Station references 


Greek 


WES 8 


gt w x 


Definition of ax/L. 

Ratio of specific heats of gas, 6% = lol. 

Change in enthalpy, (BTU/1b). 

Number of neutrons produced per absorption in the fis- 
sionable isotope. 

Definition of 054,793. 

Transport mean free path, (cm). 

Absolute viscosity of the fluid, (1b/ses ft). 

Average cosine of the scattering angle per neutron 
collision in the laboratory system. 

Number of neutrons produced per fission in the fission- 
able isotopes. 

Density, (1p/rt?). 

Microscopie cross section (em?) . 
Ta microscopic absorption eress aection (em*). 
c. microscopic scattering crosa section (em*). 


σ, microscopic fission cross section (em). 
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Macroscopic cross section (em), (ir). 
1) 


=, macroscopic absorption cross section (om 
| 1 
; 


£, macroseopie seattering eross seetion (em 
1) 


Z£, maeroscopie fission cross section ( em" 


Number of neutron absorptions / em sec. 


Tractive force at wall due to fluid friction (1b/ft?). 


Flux 1n neutrons / om^ sec. 


Laplacian operator. 








LIMITING CONDITIONS 


To obtein a realistic operating condition for the air- 
eooled nuclear heat exchanger the following conditions were 


assumed: 


M. = 2.5, h=37,500 ft, T= 392%, P« 3.09 1b/in® 


0 
Godsey and Young (1, p.221) stated that the total temperature 
through the inlet diffuser depends only on flight Mach num- 
ber and ambient temperature and is a constant regardless of 
diffuser losses if the process is considered adlabatic. 
Therefore the relation 


yal „Ê 
T0,0 * 7To,1 * 79,2 * T( * t- Mo? (1) 


may be used where the station numbers correspond to those of 
Figure 1. Godsey and Young (1, p.221) stated that the total 
pressure at the heat exchanger inlet, Po, 2* was dependent upon 
the flight Mach number, amblent pressure, free stream diffu- 
sion losses, duct friction losses, and duct diffusion losses. 
The last three effeets my be grouped under the term ram ef- 
fectiveness or pressure recovery effectiveness. From data 
presented by Godsey and Young (1, p«86), the ram effectiveness 
at M, = 2.5 was given as approximately 0.98; however, subsonie 
diffusion was assumed to be loss free. For a more realistic 


condition a value of 0.90 was chosen. The above ratios may be 


in 
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written out in the following expression for stagnation pres- 


sure at the heat exchanger inlet: 


“1 . 2, Yr-1 . Yo tual 
Pog = Pot + Sm?) o (1a) 


0,2 ideal 


Solutions to Bquations 1 and la give To,2 = BIOR and 
Po,2 * 46.5 lb per sq in. for inlet conditions to the heat 
exchanger. 

Prior to converting the above inlet stagnation pressure 
to stream pressure, it is necessary to examine the limiting 
conditions for heat exchange. Assuming a uniform cross 
sectional heat exchanger area under a given inlet stagnation 
temperature, which in this case is 86),°R, and various inlet 
Mach numbers, it 15 desired to exchange heat to the extent 
that the stagnation temperature at the exit of the heat ex- 
changor be ralsed to a maximum value. By use of methods out- 
lined in Reference (2, p.l.2l) tho results of various inlet 
Mach numbers are shown on Figure 3 which reveal that when a 
certain maximum temperature occurs, thermal choking even 
under frictionless conditions exists in the heating process 
and definitely limits the capacity of the constant area 
combustor or heat exchanger. Stated another way, if at any 
position downstream for a constant area duct, the stagnation 
temperature is equal to 70,3 maximum, the flow is chokeé 


since the Mach number is unity. Any addition of friction 


9 


will reduce the downstream interval for 70,3 maximum to occur. 
With M, = 0.2 and 70,3 maximum achieved, 5100°R, the stagna- 
tion pressure would suffer a 20 per cent loss due entirely to 
the heating effect. 

Without becoming involved with a specific shock cone and 
diffuser desien, the inlet Mach number was set at 0.2, which 


resulted in the following stream conditions: 


P, = 15.1 1b/in® 

ο» = +1385 1»/re? 

v, * 291 ft/sec 

Το x §80°R 

G = hOsk 1b/ft^ see 


As shown by Figure 2, the heat exchanger elements are 
thin hollow constant concentric area right cylinders which 
are adaptable to one-dimensional compressible gas flow treat- 
ment without introducing serious errors. This assumes that 
in any cne annulus the change of stream properties in the 
direction of flow is mich larger than in the radial direction. 
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HEAT TRANSFER ANALYSIS 


MeAdama (3, p.289) correlated the experimental data of 
Kays and London (h) and because compactness of an air-borne 
heat exchanger is an important parameter, a comparison of 
heat transfer coefficients per unit volume for a given fric- 
tion power per unit volume showed ruffled fins and Louvered 
plate fins as being the most effective. It appears that the 
highest ratio of heat transfer to pumping power is obtained 
when form drag is absent, or when all the drag is due to skin 
friction which is the case in common for flat plates or fins. 
Even though externally finned tubes are of value in increas-~- 
ing the rate of heat transfer, friction losses become pro- 
hibitive for this type of installation. For these reasons 
thin concentric hollow cylinders were chosen as the fuel ele- 
ments of the heat exchanger and were restricted to constant 


Longitudinal cross section area for nuclear simplification. 


Deve Lopment 


If heat losa and longitudinal thermal conductivity are 
neglected, the rise in the onthalpy of the air per second is 
related to the effective coefficient of heat transfer from 
an annulus by & differential eguation which is developed in 


Appendix C as 





11 
A 2 2 a a 
t (D, - Da )ev C aro = h,r(D, + D4) (To w τρ αν) dx. (2) 


The particular temperature differential is shown schematically 
in Figure lh. The temperature increment for heat transfer is 


the difference between the actual wall temperature, T » and 


DW 


the temperature T , Which the wall would assume for the 


0,aw 
same flow ias but with zero heat flux. McAdams et al. 
(5) found that the effective heat transfer coefficient, hes 
expressed in terms of the Stanton number, Noms was independent 
of the temperature potential, To,w - Tos when based on 
To,w - To,aw while the coefficients A and As were not inde- 
pendent. 

The Reynolds analogy with Nop = 1 relates the Stanton 


number to the friction factor by 


esi l (3) 


From experimental data for a moderate AT, L/D, = 60, and 
5,000 € N,.. 200,000, the following equations from McAdams 
(3, pp.155,219) may be combined 


£ = „ol6/tpe/n,)*® (4) 


by ο 2/3 . 
0 TJ. t (5) 


(26/μρ) Y 


A ο εως ὝΝ 4 
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to produce 


h 
T Af 
TT zz > (6) 


where, by definition, Az (Np) 72/3 . Since Equations 3 and 6 
are essentially the same, the limitation to a moderate tempera- 
ture change may be removed without introducing serious errors 
when using the previously defined temperature potential. The 
following relationship, where specific heat is considered a 


constant, 18 used for the remainder of the analysis: 


— (7 
p 


The friction factor of Equation 7 1s set by Equation whieh 
ía dependent upon Reynolds number and the degree of roughness. 
Since mass velocity, C, and specific heat are essentially 
constant, the friction factor may be considered to depend upon 
hydraulic diameter, viscosity, and degree of roughness. An 
average value of Af = .008 wes obtained from Equation ij for 
1/2"« Da < 1-1/2" and LhOO°R < py «1600?R. The varience of the 
friction factor within the above range of hydraulic diameter 
may be compensated for by added degrees of roughness; however, 
Cope (6) found that with water in the turbulent region even 
though friction was six times that for smooth pipes, the heat 
transfer was only 20 to 100 per cent greater than for smooth 
tubes. MeAdems (3, pe223) states that because of a very 
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small temperature change the heat transfer coefficients ob- 
tained were somewhat uncertain. Correlation with air may not 
apply due to different Prandtl number effect; however, there 
is a possibility that with small variation in hydraulic diam- 
eter the friction factor could remain a constant by varying 
the degrees of roughness. 

Experimental date have shown that for turbulent flow in 
rectangular duets and annular spaces friction may be evaluated 
from data for circular pipes using a hydraulic diameter shown 


in Figure 2 as 


a(D ©» D 2) 


D. € hr -— EI ο. - Ὁ ; (8) 
C h aD, + 1 2 l 


A combination of Equations 2 and 7 with Equation 8 results in 


— — = 2 Ar δὲ . (9) 


If a gas stream of uniform temperature is brought to rest 
isentropically, as at the stagnation of a blunt body, the tem- 


perature rise for an ideal gas will be 


2 
To - P “Ta = P(2 +f MÊ) (10) 


For convenience a temperature recovery factor Npp as shown in 
Figure lh is defined as 


πι + ne. 2 eee 
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N. s Quam I. Quam a (11) 
T 0 - T v”/2g IC, 


By a combination of Equation 10 ená 11, the stream temperature, 


T, may be eliminated to produce 


> (12) 





Substitution of Equation 12 into 9, yields 


—- ee 2Ar & o (13) 


The generalized differential Mach number equation of 
Reference (2, p.li.12) simplifica to 


(15 £m)s-L- was 
s - τὸ 


MÉ (1 + Ll) 
s re uc $E (1h) 





with negligible structure drag forces and constant specific 


heat, molecular weight, area and mass flow rate. 
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Equation 13 was substituted into Equstion 1l which gave 


ws Xu ae) 
e πῇ 


e 
T ι -κ zl μὲ 
(1 ο μ)λ] ZE.-1+ BE κος 
T 1 + Y 
0 (15) 


For a turbulent boundsry layer over a flat plate, 
MeAdams (3, p.311) recommends that the temperature recovery 
factor end the Prandtl number be related by 


N PR 


RF 
MeAdams et al. (5) reported rocovery factors for subsonic flow 
of air near the outlet and of a long smooth tube were in good 
agreement with those predicted by Equation 16. Appendix A 
shows the procedure for solving Equations 13 and 15 for one 
Mach number and stagnation temperature at an assumed position 
downstream. From these calculations it was evident that the 
contribution of the recovery factor term Nap, Wes negligible 
for Mach numbers up to approximately 0.8; however, the re- 
sulting stagnation temperatures will be nearly 1.5 per cent 
higher than those which have incorporated this term. 

For the remainder of this analysis the recovery factor 
and Prandtl number were assumed unity which reduced Equation 
13 to 
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0 ax 
— = — - 1 af x (17) 
20 0 6 


It 1s of interest to note that Equation 17 could be expected 
from experimental data by MeAdams (3, p31) because the heat 


transfer coefficient, h., deereases and approaches the constant 


g? 
value of he for temperature differences greater than 10058. 
Por a high temperature potontial, πα approaches h, and may be 
used directly in Equation 1l. | 

Equations 1l and 17 were integrated approximately over 
short intervals of Mach numbers by soda the procedure out- 
lined in Reference (2, 1.60). Appendix B contains the step- 
wise procedure used aná Tables 1-3, Appendix B, have the com- 
plete data for stagnation and stream properties as a function 
of L/D, ratios for various longitudinally constant wall tem- 
peratures. Stagnation temperatures and pressures versus L/D, 
ratios are shewn in Pigure 5. 

However, from a nuclear standpoint the longitudinal flux 
of fuel element cylindera will follow a sine distribution with 
the maximum flux produced at the midpoint. This does not 
imply that the maximum wall temperature will occur at the 
longitudinal midpoint. By use of 
= 2079 8 max - Τρ») (18) 


- 7 
0,3 0,2 1 - soc a max 
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which was developed in Aprendix C and by using the previously 
selected temperatures for maximum wall temperatures, the 
stagnation temperatures versus L/D, ratios were computed as 
shown in Table l| of Appendix C and plotted in Figure 5. Stag- 
nation pressures for variable wall temperatures were obtained 
by interpolation from previously ealeculated pressures. 

An isentrople expansion using air data from Keenan and 
Kaye (7) was assumed in order to obtain the chance of enthalpy 
from stage 3 to of Figure 1. The Brayton eycle as shown in 
Figure 6 represents an overall view; however, at various L/D, 
ratios, the change of enthalpy was computed using an ambient 
pressure of 3.09 lb per sq in. Figure 7 shows the variance 
of enthalpy for both constant and variable wall temperatures. 
The right hand ordinate of Pigure 7 related change of enthalpy 
ee sam 





4s of thrust per sq ft by 


T» a 
- = t (v - Vo) , (19) 


where Yo = 2.56 and Yh = 223.7Vh 0,3 > By, ο 
Discussion 
The calculated values of stapnation temperatures and 


pressure as shown in Figure 5 do not show the contributing 


effeets on heat transfer and friction caused by the length 


18 


to diameter ratio and the effects of sudden enlargement or 
contraction of flow. The former effect is small since the 
L/D, ratio for this heat exchanger should be greater than 60; 
the latter efrect my still be appreeieble even though the 
peripheral thickness of the cylinders are in the order of 0.05 
to 0.07 in. The specific heat ratio was set at lel, end the 
friction factor, specific heat, and mass velocity were con- 
stant throughout the heating process. The cosine flux distri- 
bution in the longitudinal direction 1s justifiable for a 
finite cylinder. 

The temeratures obteined by normal forced convection 
methods for an assumed Τον = 2)60°R are shown in Figure 5 for 
comparison. The temperatures show good correlation with stag- 
nation temperatures obtained by the Reynolds analogy but are 
considerably higher than the stream temperatures. The forced 
convection formulas will not show any decreasing trend in 
stream temperatures or a maximum stagnation temperature as 
heat is continually added as will the Reynolds method for com- 
pressible fluid which will eventually produce choked flow at 
a certain L/D, ratio downstream. 

Figure 7 shows that the effect on thrust reduction was 
quite severe for a 2000°F constant wall temperature versus a 
2000°F maximum wall temperature; the thrust per square foot 
was reduced by 18 per cent. For a uniform velocity at the 
exit of stage |, a constant change in enthalpy is required. 


If the neutron flux was such as to produce a series of maximum 
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wall temperatures from 1600 to 220.9P and with & constant 
enthalpy change of 200 BTU per 1b imposed, the UD, ratío 
would very fron L47 to 100, Sinee all cylinders sre desired 
to have the same length for nuclear considerations, the hy- 
draulic Glameters would vary by 4.7 per cent with decreasing 
hydraulic diumeters as radial distance is increased. This 
diameter variance would still be subject to the limitation 
imposed on the hydraulic diameter due to a constant friction 
factor, 

To increase the thrust which ia based on the heat ex- 
changer cross-sectional area, each fuel cylinder should oper- 
ate at the maximum enthalpy chenge ae shown by Figure 7 and 
the thrust would be a weighted average depending on the 
frontal area of the heat exchanger. With e constent length 
assumed the hydraulic diameter would vary by only 12 per cent. 
However, each annulus would be operating ata different veloc- 
ity level which may create internal vibration within the rem- 
jet pod. With a weighted radial average of 20009?P for the 
eylinder temperatures encompaseaing 7.7 sq ft of heat exchanger 
eross-sectional area, from Figure 7 the resultant thrust would 
be approximately 9,500 lb at a flight Mach number of 2.5 at 
37,500 ft. 

With inlet Mach numbers greater than 0.2, the L/D, ratio 
would be substantially reduced. From Reference (2, p.li31), 
with M, = 05, f = 04005, T, , * 10009R, P, , = 10 1b/1n*, and 
To,w = 2000°R and 150058, choked flow occurred at an L/D, ratio 
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of 20 and 35, respectively. To resist the pancaking effect 
which may reach the limit for nuclear criticality, the inlet 
Mach number appears to have a definite upper limit when heat- 
ing is desired. 


Conclusion 


Por the assumed inlet conditions, eylindrical fuel ele- 
ments and wall temperatures, the amount of heat transferred 
was substantial enough to produce 1,000 to 1,100 pounds of 
thrust per sq ft if heat exchanger frontal ares with a maximum 
wall temperatures range from 1600 to 2200°F at a flight Mach 
number of 2.5 and an altitude of 37,500 ft. With a weighted 
radial average of 2009%P for the cylinder temperatures encom 
passing 7.7 sq ft of heat exchanger area the resultant thrust 
was 9,500 lb. 
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PUEL ELEMENT CONSIDERATION 


Due to the absence of effective moderation within this 
nuclear heat exchanger, the unit may be considered as a fast 
neutron reactor and as such the restrictions normally imposed 
on structural materials by high thermal neutron absorption 
may be removed. Whether any structural material now in exist- 
ence can be used to contain, protect, and support the fission- 
able phase for the high temperatures required la the subject 
of the following development aná discusalon. 


Development 


The nuclear heat exchanger as shown in Figure 2 inecor- 
— thin concentric cylinders which sre the solid fission- 
able containers. The enriched uranium fuel cores are to re- 
main a solid to prevent serious containment problems. The 
thin sheets of the fissionable phase ere sandwiched between 
the imner and outer layers of the non-fissionable structural 
metal which predominates in volume so that the alloyed en- 
riched uranium fuel sheets are surrounded by a continuous 
network of structural metal for strength at high temperatures. 

The following are aome criteria that were considered for 


the selection of fuel elements and structural material: 
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A. Fuel element criteria 


4 
* 


5. 


Melting point above 2500?F 

Heutron's energy maintsined 

Mechanical properties simílar to the container 
material during fabrication and high tempera- 
ture operation. 

High thermal conductivity and low specific 
heat for good heat transfer 


fase of fabrication 


B. Container material criteria 


Le 


vu gf w 


High temperature mechanical properties under 
static and smell radial dynamic loads 
Resistance to oxidation at high temperatures 
High resistance to thermal shock 

Neutron energy maintained 


Ease of fabrication, rolling, and weldability. 


Discussion 


The melting point of natural uranium is 2072°F; for the 


required temperature range of 2000 to 2200°7, the unalloyed 


uranium would become a liquid and present serious containment 


problems. In order to increase the melting point, the various 


binary diagrams of uranium alloys of Saller (8, p.i13-1h) 


revealed that a uranium-zirconium system of 20.7 per cent sir- 


conium by atomic weight would have a melting point of 2h509F. 
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Elements with less than a mass number of lO were not consid- 
ered because of the increased slowing down effect on the 
neutrons. 

Zirconium alloys at high temperatures suffer a loss of 
strength and a loss of corrosion resistance. As reported by 
Dayton (9, p.l60), the strongest alloy tested at 930°F creeps 
at a rete of 10% per cent per hour with a stress of only 
4000 1b per sq in. Sinee e thin sheet is sandwiched between 
two layers of structural material, the loss of strength and 
eorrosion resistance are secondary to the requirement that the 
fissionable material remain a solid. The coefficient of 
linear thermal expansion for 20.7 per cent zirconium by weight 
was not available; however, Dayton (9, peli6l) reported a co- 
efficient for pure zirconium at 932%P for the "a" axis of 
5.8 x 1076 while that of pure uranium was anisotropic as shown 
by Saller (8, p.389). The coefficient of thermal expansion 
for this particular alloy was not found; however, from Saller 
(8, p.234) the composition for 20.7 per eent zirconium by 
weight indicated 3 phase changes for the temperature range up 
το 2200°F and whether the anisotropic effect of uranium would 
result in volume changes so severe as to make bonding or 
cladding impractical was not determined, 

The low thermal conductivities of both zirconium and 
uranium are partially compensated for by low specific heats 


and the thin sheet effect of high area to volume ratio. 
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Uranium can be hot rolled wíthout danger of cracking. 
Whether this particular alloy can be rolled into sheets with 
thicknesses as small as 0.03 inch was not determined. 

The container material for temperatures up to 2200°F 
presents a considerable problem when restricted to present day 
material availability. The wrought cobalt-base alloys, such 
as Haynes Stellite-25, Allegheny-Ludlum Steel Corporation 
3-590 and 5-816, and Westinghouse Electric Corporation Re- 
fractaloy 26 and 70, appear to have excellent corrosion re- 
sistance up to temperatures as high as 2000°F, Eichen and 
Jackson (10, p+556), in an oxidizing atmosphere. However, 
Eichen and J&ckson (10, p.539) reported the ultimate strength 
of Haynes Alloy 25 at 2200°F to be 4,000 1b per sq in. for a 
sheet 9.042 inch in thickness and annealed at 2200°F. Another 
major drawback may be the high residual actívity in the form 
of gamma radiation which ts characteristic of cobalt-base al- 
1078. 

Molybdenum is of interest as a high-temperature material 
because of its high melting point, L750?P, high strength and 
stiffness at elevated temperatures.  Northeott (11, p.106) 
reported creep-rupture strength for wrought are cast molyb- 
domm, +15 por cent titanium, of 25,000 lb per sq in at 2000?F 
for 100 hours. The thermal expansion coefficient as reported 
by Northeott (11, p.26) ranged from 5.1 x 1075 at 932%P to 
7.2 χ 107° at 3632%P, Thus the thermal expansion of both pure 


zirconium and molybdenum alloys appear to be very similar. 
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However, alloying with uranium might alter the thermal ex- 
pansion characteristic considerably. 

Uneoated molybdenum 1s useless as a material of construe- 
tion in oxidizing atmospheres above 1090°F. Norden (12, 
ppe298-30h) reported that so far no molybdenum-base alloys 
have been developed with satisfactory oxidation resistance; 
however, Hiester et al. (13) reviewed progress made with al- 
loying techniques. Sherwood (1h, p.26) stated that molybdenum 
Specimens protected by coatings 0.003 in thick and consisting 
primarily of molybdenum disilicide have succeasfully resisted 
oxidation at 3092°F for more than 1000 hours. However, the 
coating must be poreefree in order to prevent serious loss of 
material. Effects of the dísilicide coating on heat transfer 
and thermal cyeling were not included. 

The high thermal conductivity of molybdenum, .32 cal per 
see em °C, and low speeifíc hest, .061 cal per em °C, permit 
it to be heated or cooled rapidly with low resultant stress 
which give it excellent resistance to thermal shock. The heat 
transfer characteristics are excellent due to the high thermal 
conductivity. The neutron moderation effects are small since 
it has a mass number of 96. 

As brought out by Norden (12, p+298) a chief advantage 
of molybdenum over ceramics, cermets, and some super alloys is 
the variety of forms and sizes available on a production basis. 
By the arc-cast methods, sheets from 9.005 to 0.020 in. thick 
and 36 in. wide are available. Molybdenum has been joined 
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successfully by spot welding, arc welding, soldering, and 
brazing but precise controls must be used since excess heat- 


ing will result in extremely brittle welds. 


Conclusion 


Molybdenum alloy, 0.45 per cent titaniua, with a protec- 
tive disilicide coating appeared to have the best probability 
of withstanding temperatures up to 2200°F in an oxidizing at- 
mosphere besides being able to withstand thermal shock and 
offer near optimo heat-transfer characteristics. The inner 
fissionable phase of enriched uranium alloyed with 20,7 per 
cent by weicht of zirconium, even though soft and pliable at 
this temperature range, which could be an advantage as far as 
the container material is concerned, should expel the heat 
produced internally due to the high surface to volume ratio. 
Whether or not the thermal expansion coefficient of the 
urenium-zirconium would be compatible with molybdenum and that 
a metallurgical bond would be formed with molybdenum was not 
ascertained. 
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NUCLEAR ANALYSIS 


The actual heterogeneous reactor configuration as shown 
in Figure 2 was treated as a homogeneous unreflected system 
using one-group diffusion theory to obtain the critical mass 
of uranium. The flattening of flux and reduction of critical 
mass by use of a reflector were not considered. The neutron 
flux necessary to produce a maximum fuel element temperature 
was calculated. The control mechanism necessary to maintain 


the flux level was not considered. 
Development 


The diffusion equation for thermal neutrons for steady 
state conditions as developed by Glasstone and Edlund (15, 
p.192) 1s 


DIS (Fl) - 5 (F)+ pa (F, Ty) =O «. (20) 


For an idealized fast reactor composed of atoms of material 
so heavy that neutrons upen colliding rebound without losing 
an appreciable amount of energy, the slowing down density term, 
pq (r, Tip), my be considered negligible. If fast neutrons 


occur only at one energy, the source term can be written as 


3 = KES b (21) 
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Using Glasstone and Ldlund's (15) treatment the multi- 


plication factor may be simplified to 


Κον LT. 3 (22) 
=e | PUEL 
by assuming that the fast Cission factor, the resonance escape 
probability, and the fast neutron utilization are essentially 
unity. The fast neutron utilization term implies that the 
eapture of neutrons by material other than the fissionable 
phase 4s negligible. Shapiro (16, p.130) states that the ra- 
diative capture cross sections for elements having mass numbers 
less than 100 and neutrons numbers less than 60 are not greater 
than 50 millibarns, If the fast neutron utilization factor 
were not unity, the treatment given by Kaplan (17, p.525) for 
thermal neutrons would still apply. 
A combination of Equations 20 to 22 yields 


VÍA (5) + B 4 (7) =0 (23) 
where B^ s 5. (7 - 1)/0. (24) 
For a weak absorber as described by Glasstone aná Fdlund 


(15, pp.92,98), the diffusion coefficient is related to the 


transport mean free path by 
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Àt 1 
D=e-— = 3A. (25) 
where A, = A, (1 = hy) and μ 5 2/34. (26) 


As a first approximation the macroscopic transport and 
scattering cross-sections are considered to be essentially 
equal except in eases where the actual values are known, The 
atomic masses of the uranium, zirconium, and molybdenum are 
high enough to mke the Yo term of Equation 26 negligible. As 
& second approximation the elastic seattering cross section is 
expected to be equal roughly to the total cross section in the 
middle energy region of 1 to 500 kev where neutron gamna cross 
sections are small, Shapiro (16, p.70). The constants used in 
this anelysis are listed in Appendix D with additional refer- 
ences. 


Substitution of Equation 25 into 2) results in 
Bê = 38, B (N + 1) (27) 


The solution to Equation 23 for the buckling in terms of 
spatial dimensions of a right circular cylinder as developed 
by Glasstone and Edlund (15, p.213) is 


ds (2.405) a (8) ’ (28) 


where R and H are the extrapolated radius and height of the 


cylinder in centimeters, 
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A combination of Fauationa 27 and 28 ylelds 


38,E,(9 - Y = Cu + (2) (29) 
Equation 29 is the approximate critical equation for an 
idealized fast homogeneous unreflected reactor in the form of 
a right circular cylinder. Appendix DP shows the calculations 
that were used to determine the critical mass of uranium for 
a right circular cylinder with e physical diameter and length 
of h2 and 60 in, respectively. 

The various component weights which were based on volume 
fractions with corresponding densities sre calculated in Ap- 
pendix D., 

For the flux required to produce the maximum temperature 
change, the relationship 


4 (T5 To. )GC κ 
max LE St (ie) (30) 
wes used. The result is shown in Appendix D. 


Discussion 


The length of the heat exchanger was limited to 60 in. 
because Figure 7 showed that the maximum thrust per sq ft was 
obtained with an approximate L/D, ratio of 120. This produced 
a 0.50 in. hydraulic diameter or a 0.25 1n. fuel element spacing 
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The uranium-zirconium alloy with a composition of 79.3 
per cent enriched uranium by weight was selected because of 
its 2,:50°F melting point. A larger percentage of uranium 
would reduce this melting point which is considered & minimum 
for the high temperatures desirec, The alloy composition was 
considered a constant and treated as a mixture for volume 
fractions since the actual density was unknown. 

The component volume fractions of the heat exchanger 
resulted from trial and error thickness selections for the 
uranium-zirconium and molybdenum alloys which satisfied Equa- 
tion 29 for a physical diameter and length equal to or less 
than h2 and 60 in, respectively. Since the physical length 
was fixed previously for maximum thrust per sq ft, the radius 
was considered variable, 

The nuclear reactor, when treated as a homogeneous un- 
reflected reactor according to Equation 29, was determined to 
be critical at a physical diameter of lO in and a length of 
60 in. The large oritical volume, 3.6 eu ft, wes due to the 


low macroscopic absorption of uranium, 0.0056 em > 


; Which re- 
sulted indirectly from the 75 per cent volume of air. A dè- 
crease in the volume fraction of uranium would decrease both 
the total macroscopic absorption and transport eross sections 
with a resultant increase in radius if the length remained 
constant. 

The weight of enriched, 90 per cent uranium required for 


the critical volume of 43.6 cu ft was 3,622 lbs. The amount 
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of zirconium and molybdenum was 9609 and 3,1470 lbs, respective- 
ly, which resulted in a total heat exchanger weight of 8,052 
lbs. At current prices of approximately 30 dollars per gram, 
50 million dollars would be required for the 90 per cent en- 
riched uranium. 

The flux that was necessary to generate a maximum wall 
temperature of 2200°F with a 0.0315 inch fuel element thick- 
ness was 6 x 1014 n per sq em sec which assumed zero radial 
flux depression across the fuel element. 

Since the radial flux will follow an approximate cosine 
distribution, 1t is apparent that for an average maximum wall 
temperature of 2000°F which was assumed to obtain 9,500 lbs 
ef thrust with a 7.7 sq ft frontal area, the fuel elements 
near the centerline would operate at temperatures exceeding 
2200°F. A proposed central control rod and reflector would 
tend to flatten the flux; however, the flux envelope resulting 
from these effects was not investigated. 

The application of the one-energy group diffusíon theory 
to this analysis was a realistic approximation because the 
transport mean free path, 8.7 em, was small in comparison to 
the boundary distence. The extrapolation distance, Rer, 48 
not generally equal to 0.71 A, for a fast reactor and its 
exact value can be obtained only by complicated methods, Glas- 
stone (20, p.217). However, as pointed out by Classtone (290, 
p+217), there is probably no point in a fast reactor that is 


more than two scattering or transport mean free paths from a 
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boundary. This statement 1s supported by the large quantities 
of fuel required for the heat exchanger which has a minimum of 
11 transport mean free paths in the radial direction. 

An intermediate neutron energy reactor may offer the pos- 
sibility of reducing the fuel quantity and weight while still 
maintaining a large volume fraction of air for compact heat 
transfer. The slowing down density term of Equation 20 was 
neglected because air, which composed 75 per cent of the heat 
exchanger volume, has a negligible moderation effect on neu- 
trons. To continue this simplification the material for the 
fuel alloy and container was selected from elements with rela- 
tively high mass numbers. However, with the proper selection 
of lighter materials for the fuel alloy, container and reflec- 
tor, the moderation effects may increase to the extent that 
this fast reactor would evolve into an intermediate reactor. 
It is believed that any trend toward a thermal reactor would 


decrease the fuel requirements for the configuration. 


Conclusion 


From a nuclear aspect this heat exchanger, which was com- 
posed by volume of 12.5 per eent alloyed fuel, 12.5 per cent 
molybdenum structure, and 75.0 per cent air, was critical at 
a diameter and length of 0 and 60 inches, respectively. The 
total weight was 8,052 lbs which included 3,622 lbs of enriched 
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90 per cent uranium. Approximately 50 million dollars would 
be required for this amount of uranium. 

The flux that was necessary to generate a maximum wall 
temperature of 2200°F with a fuel element thickness, 0.0315 
in. was reasonable at 10!?n per sg cm sec, 

With cost as a major consideration there 1s probably no 
point in a fast reactor more than several scattering or trans- 
port mean free paths from a boundary. An intermediate neutron 
energy reactor may offer the possibility of reducing the fuel 
guantity and weight for this configuration. 





OVER-ALL PEASIBILITY CONCLUSIONS 


For an air inlet velocity of M = 0,2, & weighted radial 
average of 2000°F for the maximum cylinder temperature, a 
frontal area of 7.7 sq ft, a h2 in diameter, and a 60 in 
length, the amount of heat transferred was substantial enough 
to produce 9,500 lbs of thrust with a flight Mach number of 
245 at 37,500 ft. 

In order to obtain an average temperature of 2000°F the 
fissionable material within the thin right circular cylinders 
hed to reach a critical condition. This required 3,622 lbs 
of enriched 90 per cent uranium with an approximate cost of 
50 million dollars. The heterogeneous nuclear heat exchanger 
was critical with a diameter and length of ΠΟ and 60 in, 
respectively and a total weight of 8,052 lbs. 

The low thrust to weight ratio may be compensated for by 
the negligible fuel burn up when compared to a turbojet sys- 
tom. 

Molybdenum alloy, 0.45 per cent titanium, which was pro- 
tected with a disilicide coating appeared to have the high 
temperature properties for the container material. Sandwiched 
between the molybdenum was an alloy of uraniumezirconium which 
contained 79.3 per cent uranium by weight. This binary system 
was selected because of the 2):50°F melting point at the above 
composition. Whether a metallurgical bond could be formed be- 
tween the molybdenum and fuel alloy was undetermined, 
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From the above theoretical conclusions and limiting 
conditions, the over-all feasibility of a compact axial flow 
nuclear heat exchanger was favorable for thrust, metal avail- 
ability, and nuclear criticality, borderline for thrust to 
weight ratio, end unfavorable for the cost of fuel. However, 
the proposed central control core and peripheral reflector 
should inerease thrust and reduce the cost of fuel because of 
savings in critical mags. 

An intermediates neutron energy reactor may offer the pos- 
sibility of reducing the fuel requirements. With the proper 
selection of lighter materials for the fuel alloy, container, 
and reflector, the moderation effects may increase to the ex- 
tent that this fast reactor would evolve into an intermediate 


reactor. 
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Figure 2. Nuclear heat exchanger. 
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Relationship of various temperatures end heat 
transfer coefficients where T, » T = V^/2gJC 
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Figure 6. The Brayton cycle for temperature and entropy 


conditions for stations corresponding to Figure 1, 
Values for dry air are from Reference (7). 
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APPENDIX A 


Stepwise Integration for Stagnation Temperature 


end Mach Number as a Function of L/D, Ratio 


The procedure used for integration of Equations 13 and 
14 are outlined in Reference (2) and are reproduced in part 
using the following approximations: 
1. The integral of aT To may be approximated by 
2(7, o - 70,17 Toa * Το»). 
2. The coefficients of dx/D, in Equations 13 and lih are 


taken at M2 = (M, + Μο) /8. 


The inlet conditions are as follows: 


My = 0.2 Nop = 0.66) 

To,1 = BEL®R ap = 0.872 

Po * hó.5 1b/1n* A s 1.32 
> 

Poy * 2h60*R f = 0,008 


For a section 2 at a distance o “E, downstream from section 
1, the following computations are carried out. 

a. Let (x, - x,)/D, = 30 

b. Estimate M, = 0+3, Μι ϱ = 0.25 


δ. Calculate ao 2 > Τρ, 1) Too " 70,1) from Équation 


13. 





ly? 


0,2 ^ To,1! To (1 = Hop) 2 
2 p—LU-* | Tit Ha Y | 
0,2 0,1 1 - > M 2 
ϑ 


0,2 + 0,1 1 


2AÀf 2 1 


9 


(To ΙΙ 88h) 60 1 - ¿872 2) (0625 
- “α-.--ᾱ, 7 - Br 1 + .2(.062 F 
0,2 + 88h 0,2 εδί. 
2(1,32) (.008) (30) 
70,2 * 16h0 


d. Caleulate M, from Equation 15 and compare with as- 
sumed value of My 


3 . 06 + 


M - M^ x (1 = 0685] 06 2(4008) (30)* 


fa + 1.4 x 0625) (1.32) | .95) *2x Lok sonas] 


M, a Μις = 0.0500 


n, = 0,0500 + 0.0100 


N, = 0.3 This agrees with the assumed value, 





e. If atep å was not : 
atep s with an im wales of T 


al, 

f. After convergence is obtained, the process 1s re- 
"c — ow ot - mra s 

| peated to find the Mach number at a section 3. 

m" ο” ο... 4 Y Da 
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APPENDIX B 


Stepwise Integration for Stagnation Temperature 
and L/D, Ratio as a Function of Mach Number 


The procedures used for integration of Fquations 1 and 
17 are outlined in Reference (2). Eaustion lh may be inte- 
grated aprroximately over a short interval of Mach number by 
assuming that the influence coefficients have a constant value 
corresponding to (M, + M,)/2, that the magnitude of T, is 
eonstant during integration with a value of To 1 + 70,2)/?: 
and that the frietion factor 1s constant. With these approxi- 
mations and by the substitution of Equation 17 into 1h, Eque- 
tion 1h integrates into 


T δν ar 


u: (17a) 


For the same inlet conditions as were given in Appendix A, 


let 
ELE x, = 0.0576 
M, = 0420 μις = 0,000 


TA S | M, = M? = 0,0176 





50 


From Reference (2), P, = 0.05483 and P, = 0.00348. Heuation 
ο 


17a converged after three attempts to 


T 


0 | 
vas * 1.3hh or T, , * 1188*R 


, 


Equation 1l gave 


qam la 
X, = X 
2 1 0,1 1 785 - 
m p m py = y "^ "- = 16.85 
e ar 04H _ 0,2 n " , 
0,1 0,1 


Since both Mach number and stagnation temperatures are known 
at each section, the following equations were used for comput- 


ing the remaining properties of the air: 
- 2 
T T (1 + Ja ) 
- 0,2 
" | = 1.34 
Y= 
1 0,1 (1 € M, ) 


νὰ - 


- = 0.965 


CE 


< 
je 


< 
pt 
pwd 
Ah? 


= 1.392 


7 


Mies © om NR e vun men em 












E 1 A 
— Eze nd ada 


bacs "m MS rz 





- E 
v BEEN ΜΝ 
ι 
T Fw 
— -- D. 
9,» 
Ha > l | 
4 













—— 0 








e 
> ps 


— 


= 
κε 
- ug 





51 


£2. 
ei 


-1 „2 y 
Poe.’ TR 71 
0,1 1 1 + 5 My 
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Tables 1-3 which fellow have calculate values for the various 
properties of air as a function of L/D, ratios and different 


constant wall temperatures. 


Table 1. Air properties as a function of Mach number for a 
nF ne of 2660°R and a friction 
M To νο. τ P V e Py 
en Ἢ 1b/in^  ft/see  1b/ft?  1b/1n* 
0.20 88), 0 68ο 45.1 291 0,1385 h6.5 
0.24 1192 11.90 1182 13.6 kos 0,1000 45.0 
0,28  1h6h 25.50 1h66 hl. 5 524 0.0770 43.0 
0,32 1730 9.72 1709 39.2 646 0.0625 10.6 
0.36 1939 56.04 1892 36.8 766 0.0528 38.1 
010 2083 69.94 2020 35.3 881 0. 0461 35.6 
Ook 2182 81.74 2104} 31.8 990 0.0410 33.0 
0.48 2252 91.55 2158 29.6 1097 0,0372 | 31,8 Ἢ 
0,52 2299 99.19 2102 27.3 1188 0.0340 30.4 
0,60 2356 109.52 2182 23.6 1370 0.029) 28.6 








Table 2. 


0,20 
0.24 
0.28 
0,32 
0.36 
0.110 
0.144 
0.48 
0.52 
0,60 
0.70 
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Air properties as a function of Mach number for a 


constant wall temperature oi 21609?R and a friction 


factor of .908 

To ο, T 

on 95 
88h 0 880 
1184  13.h3 1180 
1h73 29.32 1h62 
1711 46.70 1692 
1891 63.98 1863 
2018 79.33 1970 
2100 92,25 2050 
2155 102.6 2080 
2192 110.7 2097 
2236 121.9 2100 
2258 128.2 2050 


P 
1b/in* 


45.1 
43.4 
h1.5 
38.9 
36.1 
33.7 
31.2 
29.9 
26.8 
23.2 
19.9 


V 
ft/seo 


291 
hos 
527 
649 
768 
078 
985 
1085 
1182 
1363 
1590 


S 


lb/ft? 


0,1385 
0.0998 
0.0769 
0. 0625 
0.0533 
0. 0465 
0. 0515 
0.0376 
0.0346 
0.0298 
0. 0259 


Po 


10/ in? 


16.5 
ll. 9 
h2.9 
40.3 
37.6 
34.8 
32.8 
3164 
30.2 
27.8 
25.8 








Table 3. 


.20 
n 
.28 
.32 
» 36 
oO 


46 
052 
«60 
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Air properties as a function of Mach number for a 
constant wall temperature of 2060°R and a friction 


factor of .008 
To L/D, T P Ρο 
OR oR 1b/in* 1b/1n* 
884 0 880 145.1 46.5 
1178 17.95 1160 ¡CAN lli. 6 
15437 39.75 1h11 41.0 12.2 
1631 63.55 1592 38.2 39.4 
1760 85.95 1711 35.2 36.3 
1839 104.7 1775 32.3 33.3 
1886 102,0 1811 29.6 31.1 
1911 131.9 1830 26.9 28.8 
1935 140.6 1530 24.3 26.9 
19145 145.6 1800 21.1 25.3 
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APPENDIX C 


Calculation for Outlet Stagnation Temperatures as a 


Funetion of L/D. Ratios and Maximum wall Temperatures 


with reference to Figure 2, for a heat balance on the 


fwel cylinders without axial conduction requires that 


. (4) 





The longitudinal volumetric heat source term 18 zero when 

x= 0, or L, which assumes zero extrapolation distance. 
Neglecting axial conduction within the fluid and any heat 

losses, a heat balance on the air in the annuli requires that 


pvc D aT 


Tou ~ To * —— m ' (B) 


Equating Equations A and B and performing the required 
integration results in 


5-152 * ge (st) | 1 - cong |. (c) 


Adding Equations A and C to eliminate the stagnation tem- 
perature of the air yields 











max πχ 
Tow ^ *o,2 " 7 2h ντ 
2G t 
max ( L ) | i E 
x G C De L 


in 


"s 


e) an] 
(5) 


From Equatíon E there appears to be a maximum which repre- 
sents the highest temperature of the fuel cylinder in the 
longitudinal direction. The point at which the maximum is at- 
tained is found by differentiating Equation E and setting the 


result equal to zero, 


D 
Xx * = are tan (- 5 2) (P) 


The maximum wall temperature occurring in the fuel eylin- 
der can now be obtained by inserting Equation ? into equation 
E. For purposes of computation, it is convenient to represent 


ax/L by the symbol a. 


a E xx/L (a) 
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By substituting G into Equation P, results in 


D 
tan a. * - Στ T (H) 


c. 


Utilizing Equation H, Equation E reduces to 


To,w " "es * re ($) | © Ἡ (1) 


’ max p 
From Equation 17, which was previously derived, 


aT 


ae =e 
Ow- ο e 


and introducing Equation A, which was símplified by substitut- 
ing h = roc /e, the following differential equation was de- 


rived: 
at, = gap” sin dx (3) 
Integrating Equation J results in 


maz e Γι 
To,3 Το," "E: (>) (K) 


However, Equation K may be combined with Equation I to elimi- 
nate Q naX t, wbich results in 
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al di 11... 
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s rá 2 ΟΜ αχ 0,2 
0,3 0,2 1 - BOCH, ax (18) 


The magnitude of To utilized the previously se- 
"max 
lected wall temperatures. Since T, „ was the given inlet 
8 
conditions, a4, was computed for various ιο ratios result- 
ing ín 70,3 as a function of Tom and aa * Teble 1 of 
Appendix C presents caleulated values of To 3 which are shown 
2 


in Figure 5. 


Table lh. Outlet stagnation temperatures as a function L/D, 
ratios and maximum wall temperatures 








L/D, ho 60 80 100 125 140 

tena ax -11.910 -3.270 -2.459 -] «965 -1.635 -l.e 405 
sect, ax «ὃν OLO -3.420 -2e 650 “22219 “1,920 -1./28 
l-seea,. X 6,010 h. 20 3,650 3.210 2.920 2.728 


ιο 
79,2 = 88h°R and Tax = 2660°R 











79,3 oR 1h75 1688 1858 1990 2101 210% 
T) 2 = BEKOR and TO = 2h609R 
70,3 °R 10686 1597 1748 1865 196h 2038 
= ιο - 20609 
T0,2 BEHOR and S ar 20609R 
T 9R 1275 1h15 1529 1616 1689 1715 


0,3 
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APPENDIX D 
Critical Calculations for a Homogeneous Unreflected Reactor 


The effective neutron energy was taken as 0.3 Mev based 
on flux spectra computed for existing fast reactors by Murray 
(18, pe69). From the same Reference the following constants 


were used: 
υ238 ι 7, = 0.18 barns — C, * 0.0 7, = Tel barns 
υ235 ' Ta = 1.78 barns To = 1,50 barns σ, - T», barns 


Ve 2.51 


Prom Hughes and Harvey (19) at 0.3 Mev, the total cross 


sections of the selected reactor materiala and aír were as 


follows: 
Mo 1 P 9,0 barns 
Wo: o = y. 0 barns 
: CO. = 3.0 barns 


+3 


The above total cross sections were considered to be equal ap- 
proximately to the scattering or transport cross sections. 
A 90 per cent enriched uranium mixture by weight produced 


235 N. 


= .90 or "8 = 0.105 . 
238 Mg + 235 Ne x 
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The solution to Equation 22 with the above result ylelded 


σι 


7 = = 1.50 
+ 
C'a,5 ^ Fast (Ν/Α) 


«αν... = 
«5 1,78 + (.18) (105) 


2.09 . 


Since the alloyed density of the uranlum-zirconium system 
was unknown, the alloy, 79.5 per cent uranium by weight, was 


treated as a mixture to obtain the following volume fractions: 


for = 0.432 


A hyéraulie diameter, D, = 0.5 in, was obtained with a 
selected L/D, = 120 and L = 60 in. The following thicknesses 
were assumed: 

Urenium-zirconium alloy: 0.93125 in 
Double wall of molybdenum: 0.03125 in 
Hydraulic diameter/2: 0.25000 in 
Thickness of air spaces 0.18750 in 


The approximate component volume fractions obtained from the 


above solidity factors were as follows: 


0.5925 
0,1575 
f., = 0.05h0 

Mo 0.1250 
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The macroscopic cross sections were computed for each of 
the above volume fractions which gave for a total 


E, * 0,1183 m and — £, 0.0086 . 


Solution to FTeuation 27 yielded 
B^ » 21,7 x 1074 me, 


With a physical length equal to 60 in or 152.) em and the 
transport mean free path equal to 6,15 em, the length includ- 
ing the extrapolation distance was 


H = h + TIA, = 15204 + 6.0 = 15001, om > 


The solution to Equation 29 for the critical radius in- 
cluding the extrapolation distance was 56.7 em; the physical 
radius was 50.7 em or 20 in. 

Since the total volume was 12.35 x 10? eu em, the follow- 
ing component weights were determined by using the previously 
computed volume fractions with corresponding densities: 


σ338 362 lbs 
233 3260 lbs 
Zr 960 lbs 


Mo 3470 lbs 


Total Weight 8052 lbs 
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From Equation K, Appendix C, the flux that was necessary 
to generate a maximum wall temperature of 2200°F with a fuel 


element thickness, 0.0315 in, was obtained with 


μας 7 | | | 
μμ δ. Ὁ (1/D,) 


zi 
it 


200 Mev/fission 
f 0,045 em 
1652 x 10° 
ho. 1b/ft^ sec 
.2671 BTU/1b °F 

L/D, = 120 
- To 2 = 2101 = 884° fron Figure 5. 


M 
ii 


16 ary 


Pal 
4 Y o 4 


Bo, 3 


The flux obtained from the above equation was 6.25 x 10H / 


em? sec. 
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